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Abstract 

Effects of supersymmetric particles on flavor changing neutral current and 
lepton flavor violating processes are studied in the supersymmetric SU(5) 
grand unified theory with right-handed neutrino supermultiplets. Using in- 
put parameters motivated by neutrino oscillation, it is shown that the time- 
dependent CP asymmetry of radiative B decay can be as large as 30% when 
the r — > fi 7 branching ratio becomes close to the present experimental upper 
bound. We also show that the B s —B s mixing can be significantly different 
from the presently allowed range in the standard model. 
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Although the standard model (SM) of the elementary particle theory describes current 
experimental results very well, particles and interactions outside of the SM may appear be- 
yond the energy scale available at current collider experiments. One of indications is already 
given by the atmospheric and the solar neutrino anomalies which have been interpreted as 
evidences of neutrino oscillation A natural way to introduce small neutrino masses 

for the neutrino oscillation is the see-saw mechanism || where the right-handed neutrino is 
introduced with a very heavy mass. This scenario suggests the existence of a new source of 
flavor mixings in the lepton sector at much higher energy scale than the electroweak scale. 

In this letter we consider flavor changing neutral current (FCNC) processes and lepton 
flavor violation (LFV) of charged lepton decays in the model of a SU(5) supersymmetric 
(SUSY) grand unified theory (GUT) which incorporates the see-saw mechanism for the 
neutrino mass generation. In the SUSY theory, the superpartners of quarks and leptons, 
namely squarks and sleptons respectively, have new flavor mixings in their mass matrices. In 
the model based on the minimal supergravity these mass matrices are assumed to be flavor- 
blind at the Planck scale. However renormalization effects due to Yukawa coupling constants 
can induce flavor mixing in the squark/slepton mass matrices. In the present model sources 
of the flavor mixing are Yukawa coupling constant matrices for quarks and leptons as well 
as that for the right-handed neutrinos. Because the quark and lepton sectors are related by 
GUT interactions, the flavor mixing relevant to the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix can generate the LFV such as fi — > and r — > /X7 processes [§J in addition to 
FCNC in hadronic observables ||. In the SUSY model with right-handed neutrinos, it is 
possible that the branching ratios of the LFV processes become large enough to be measured 
in near- future experiments . When we consider the right-handed neutrinos in the context 
of GUT, the flavor mixing related to the neutrino oscillation can be a source of the flavor 
mixing in the squark sector. We show that due to the large mixing of the second and 
third generations suggested by the atmospheric neutrino anomaly, B s -B s mixing, the time- 
dependent CP asymmetry of the B — > M s 7 process, where M s is a CP eigenstate including 
the strange quark, can have a large deviation from the SM prediction. 
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The Yukawa coupling part and the Majorana mass term of the superpotential for the 
SU(5) SUSY GUT with right-handed neutrino supermultiplets is given by 

W = y^l'M'jH:, + jy^^jih, + fSN&Hs + ^M^NiNj , (1) 

where and iVj are 10, 5 and 1 representations of SU(5) gauge group. i,j = 1, 2, 3 are 

the generation indices. H 5 and Hi are Higgs superfields with 5 and 5 representations. In 
terms of the SU(3) xSU(2) L xU(l)y, contains Q< (3, 2, |), U t (3, 1, -§) and E { (1, 1, 1) 
superfields. Here the representations for SU(3) and SU(2) groups and the U(l)y charge 
are indicated in the parentheses. includes Di (3, 1, |) and Lj (1, 2, — |), and iVj is a 
singlet of SU(3)xSU(2)ixU(l)y. M u is the Majorana mass matrix. Below the GUT scale 
(si2x 10 16 GeV) and the Majorana mass scale (= Mr) the superpotential for the minimal 
super symmetric standard model (MSSM) fields is given by 

Wmssm = JVQ,l',H> + fSQiDjH, + ./•;;'/•:,/.,//, + hH x H 2 - ^<(L^ 2 )(L i i? 2 ) , (2) 

where k v is obtained by integrating out the heavy right-handed neutrino fields. At the 
right-handed neutrino mass scale k u is given as n l J = (/^M" 1 /^)*-'. The Yukawa coupling 
constants f{j, f% and fl are related to the coupling constants ftf and fp at the GUT 
scale. The quark, charged lepton and neutrino masses and mixings are determined from the 
superpotential Eq. @ at the low energy scale. 

As discussed above, the renormalization effects due to the Yukawa coupling constants 
induce various FCNC and LFV effects from the mismatch between the quark/lepton and 
squark/slepton diagonalization matrices. In particular the large top Yukawa coupling con- 
stant is responsible for the renormalization of the qi and ur mass matrices. At the same 
time the mass matrix receives sizable corrections between the Planck and the GUT scales 
and various LFV processes are induced. In a similar way, if the neutrino Yukawa coupling 
constant is large enough, the II mass matrix and the cLr mass matrix receive sizable 
flavor changing effects due to renormalization between the Planck and the Mr scales and 
the Planck and the GUT scales, respectively. These are sources of extra contributions to 
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LFV processes and various FCNC processes such as b — > 57, the B°-B° mixing and the 
K°-K° mixing. 

It is particularly interesting that the chiral structure of the FCNC amplitudes due to 
the dpi mixing is different from that expected in the SM. For example, the flavor mixing in 
the d,R sector generates a sizable contribution to the b — > s 7^ amplitude through gluino-c/ij 
loop diagrams, whereas this amplitude is suppressed by a factor m s /rrib over the dominant 
b — > s7l amplitude in the SM. When the amplitudes with both chiralities exist, the mixing- 
induced time-dependent CP asymmetry in the B — > M s 7 process can be induced. Using 
the Wilson coefficients cj and c' 7 in the effective Lagrangian for the b — > 37 decay C = 
c-jsa^buFav + c'jSa^biF^ + H.c, the asymmetry is written as 



where T(t) (T(t)) is the decay width of B°(t) -> M s 7 {B°(t) -> M.7) and M s is some CP 
eigenstate (£ = +1(— 1) for a CP even (odd) state) such as -K^tt 0. Am d = 2\Mi 2 (B d )\ and 
#b — ar g^i2(-Bd) where M 12 (-Bd) is the B^-B^ mixing amplitude. Because the asymmetry 
can be only a few percent in the SM, a sizable asymmetry is a clear signal of new physics 
beyond the SM. 

We calculated the following observables in the FCNC and LFV processes: the CP vi- 
olation parameter in the K°-K° mixing Ek, B^-Bcl and B s -B s mass splittings Am^ and 
Am s , respectively, A t and the branching ratios B(6 — > 57), B(/x — > ej), B(r — > 7x7) and 
B(r -> e7). We solved renormalization group equations (RGEs) for Yukawa coupling con- 
stants and the SUSY breaking parameters numerically keeping all flavor matrices. After 
demanding the condition of radiative electroweak symmetry breaking, the free parameters 
in the minimal supergravity model are the universal scalar mass m , the universal gaugino 
mass Mq, the scalar trilinear parameter A , the ratio of two vacuum expectation values 
tan p and the sign of the Higgsino mass parameter /j. In addition we need to specify 
neutrino parameters. The phenomenological inputs from neutrino oscillation are two mass- 
squared differences and the Maki-Nakagawa-Sakata (MNS) matrix. In order to relate these 
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parameters to and M v , we work in the basis for Ni, Li and E\ where f l L 3 = f L J and 
In = (/atVl) 1,3 (/l an d /at are diagonal matrix) at the matching scale Mr. In this basis 



Hi, 



/nM^ 1 /nVl = ^mns^^^mns where ^m NS is the MNS matrix at Mr and k v is a 



diagonal matrix. Note that although Vl = V^ NS when M v is diagonal in this basis, two 
are independent in a general case. Once we specify three neutrino masses, Vmns, Vl and 
/jv we can determine the M v matrix. Then using the GUT relation for Yukawa coupling 
constants, we can calculate all squark and slepton mass matrices through RGEs. Note that 
Vl essentially determines the flavor mixing in the d R and II sectors in this basis. 

As typical examples of the neutrino parameters, we consider the following parameter sets 
corresponding to (i) the Mikheyev-Smirnov-Wolfenstein (MSW) small mixing angle and (ii) 
the MSW large mixing angle solutions for the solar neutrino problem ||. 
(i) small mixing: 



m. 



V 



MNS 



2.236 x 10~ 3 , 3.16 x 1(T 3 , 5.92 x 1(T 2 J eV , 
/ 0.999 0.0371 X 



-0.0262 0.707 0.707 
0.0262 -0.707 0.707 



(ii) large mixing: 



m. 



(4.0 x 10~ 3 , 5.831 x 10" 3 , 5.945 x 10~ 2 ) eV , 

( 1/V2 1/V2 ^ 



V 



MNS 



V 



/ 



-1/2 1/2 1/V2 
1/2 -1/2 l/y/2 

In each example we also take M v to be proportional to a unit matrix with a diagonal 

mns an d In 



element of M R = 4 x 10 14 GeV so that V L = K ot 



M R k%. We fix m 



pole 



175 



GeV, m 



pole 



4.8 GeV and the CKM parameters as V cb = 0.04, \V ub /V cb \ = 0.08 and take 



several values of the phase parameter in the CKM matrix 5 13 ||. We take tan/3 = 5 and 
vary other SUSY parameters mo, Mo, Aq and the sign of \x. Various phenomenological 
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constraints from SUSY particles search are included (for detail see ||). We also impose 
2 x 1CT 4 < B(6 ^57) < 4.5 x 1CT 4 |Ll| in the following analysis. 

Let us first discuss B(/i — ► ej) and Ek which turn out to be strong constraints on 
the parameter space in this model. Fig. [1] shows the correlation between B(/z — > ej) and 
B(r — > fi'y) for the neutrino parameter set (i) and (ii) for 5 13 = 60°. We can see that the 
B(/i — > ej) becomes a very strong constraint for the case (ii), which is a reflection of the 
large 1-2 mixing in the Vmns matrix. By requiring B(jU — > ej) < 1.2x 10~ n [pT2l , B(r — > ^7) 
becomes less than 10 -8 for the case (ii) whereas it can be close to the present experimental 
bound (1.1 x 10~ 6 |13|) for the case (i). We also calculated B(r — > cy) which turns out to 
be smaller than 3 x 10~ 12 in both cases. The constraint from Ek depends on the parameter 
5i 3 . After imposing the B(/i — > ej) constraint, Ek can be enhanced by 50% for the case (i) 
and by a factor of 2 for the case (ii). This means that compared to favorable values in the 
SM (50° < 5i3 < 90°), a smaller value of £13 is allowed due to the extra contributions. 

The upper part of Fig. ^| shows a correlation between Am s / Am^ and B(r — > fij) for 
case (i) and #13 = 60°. Here we imposed the constraints from B(/i — > cy) and Ek- We 
also imposed the constraint from Am^ itself though the deviation of this quantity from the 
SM value is within 5%. For the theoretical uncertainties we allow ±15% difference for Ek 
and ±40% for Am^. For Am s / Am^ we fix the hadronic parameters as /s s //s d = 1.17 and 
B Ba /B Bd = 1. We can see that Am s /Am d can be enhanced up to 30% compared to the 
SM prediction. This feature is quite different from the minimal supergravity model without 
the GUT and right-handed neutrino interactions |5) where Am s / Arria is almost the same 
as the SM value. A t for the same parameter set is shown as a function of B(r — > fi'y) in 
the lower part of Fig. |2|. We can see that \A t \ can be close to 30% when B(r — > fij) is 
larger than 10~ 8 . The large asymmetry arises because the renormalization effect due to 
/jv induces sizable contribution to c' 7 through gluino-rf fi loop diagrams. The corresponding 
figure to Fig. ^ for the case (ii) shows that B(r — » ^7) is cut off below 10~ 8 and the maximal 
deviation of Am s / Awid from SM is within 6% and \A t \ becomes at most 6%. 

In Fig. [| we show Am s / Arrid for several values of #13 for the case (i) and (ii). In these 



figures we impose B(/x — > cy) constraint. Thin vertical lines correspond to the case without 



the experimental constraints from Er, Am^ and the lower bound for Am s /Am,j ||10|| , and 
the thick lines are allowed range with these constraints. The allowed range of the SM is 
also shown in these figures. Because the new contributions to B^-B^ amplitude is small, 
the time-dependent CP asymmetry of B — > J/ip K$ in this model is essentially the same as 
the SM value, therefore we can obtain information on 5^ once this asymmetry is measured 
experimentally. For example the asymmetry of B — > J/ipKs mode is 0.4 for 813 = 25° and 
0.65 for 813 = 75°. This figure means the possible deviation from the SM may be seen in 
both cases once the CP asymmetry of B — > J/tpKs mode and Am s /Amd are measured. 

In our example we took M R which corresponds to the upper bound of fx because a 
larger Mr would lead to the blow-up of below the Planck scale. If we take a lower value 
of Mr the flavor changing amplitudes essentially scale as Mr. 

Finally we would like to comment on genarizations of our result. Firstly for a large 
tan/? case, the constraint from B(/x — * cy) becomes stronger. As a result we cannot see 
any deviation from the SM for the large mixing case in the figures corresponding to Fig. 
and |3] for tan (3 = 30, whereas the result is similar for the small mixing case. Secondly 
we considered the vacuum oscillation case. We see that the pattern of the deviation from 
the SM is similar to the small mixing case because the effect of the flavor mixing in 1-2 
generations turns out to be suppressed by the degeneracy of the two light neutrino masses. 
In this letter we took the case where M v is proportional to a unit matrix. When we consider 
a more general case, Vl is not necessarily equal to V^ns- m addition, the superpotential flU) 
does not lead to the realistic fermion mass relation especially for the first and the second 
generations. In order to solve this problem we may have to introduce more free parameters. 
Although the precise values of the predictions depend on the detail of the model in question, 
the possible new physics signals such as B(r — > 7x7), B s -B s mixing and A t may be expected 
as long as some of the neutrino Yukawa coupling constants are large. Because these signals 
provide quite different signatures compared to the SM and the minimal supergravity model 
without GUT and right-handed neutrino interactions, future experiments in B physics and 



LFV can provide us important clues on the interactions at very high energy scale. 
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FIG. 1. The correlation between B(/i —* cj) and B(r — > ^7) for small and large mixing angle 



cases. See text for the input parameters. 
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FIG. 2. The ratio of B s —B s and B^-Bd mass splittings Am s / Arrid and the magnitude factor 
At of the time-dependent CP asymmetry in the B — > M s 7 process as a function of B(t — > fi 7) for 
the small mixing case (i). 
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FIG. 3. Possible range of Am s / Arrid as a function of 5i3 for small and large mixing cases. Each 
thick vertical lines shows allowed range with the experimental constraints (see text). Thin vertical 
lines correspond to the case without the constraints from ejc, Arrid and Am s / Am^. Horizontal 
line shows the experimental lower bound of Am s / Am^- The dotted line corresponds to the SM 
value and the solid section shows the allowed range in the SM. 
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